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Abstract 

A dust cloud of Ganymede has been detected by in-situ measurements with the dust 
detector onboard the Galileo spacecraft. The dust grains have been sensed at altitudes 
below five Ganymede radii (Ganymede radius = 2,635 km). Our analysis identifies the 
particles in the dust cloud surrounding Ganymede by their impact direction, impact 
velocity, and mass distribution and implies that they have been kicked up by hyperve- 
locity impacts of micrometeoroids onto the satellite's surface. We calculate the radial 
density profile of the particles ejected from the satellite by interplanetary dust grains. 
We assume the yields, mass and velocity distributions of the ejecta obtained from lab- 
oratory impact experiments onto icy targets and consider the dynamics of the ejected 
grains in ballistic and escaping trajectories near Ganymede. The spatial dust density 
profile calculated with interplanetary particles as impactors is consistent with the pro- 
file derived from the Galileo measurements. The contribution of interstellar grains as 
projectiles is negligible. Dust measurements in the vicinities of satellites by spacecraft 
detectors are suggested as a beneficial tool to obtain more knowledge about the satellite 
surfaces, as well as dusty planetary rings maintained by satellites through the impact 
ejecta mechanism. 

Keywords: a. Asteroids, comets, meteoroids 

1 Introduction 

Since December 1995 the Galileo spacecraft has been on a bound orbit about Jupiter and 
performs regular close flybys of the Galilean satellites. With the dust detector onboard 
(Griin et al., 1992), at least three populations of dust have been identified in the Jo- 
vian system (Griin et al, 1997, 1998): 1) Streams of ten-nanometre grains were detected 
throughout the Jovian magnetosphere and were recognizable even in interplanetary space 
out to 2 AU from Jupiter (Griin et al, 1993). Recently, it has been shown that these 
streams originate from Io, the ultimate source probably being the most powerful of its 
volcanic plumes (Graps et al., 2000). 2) Bigger micrometre-sized particles form a tenuous 
dust ring between the Galilean satellites (Colwell et al., 1998; Thiessenhusen et al., 2000; 
Krivov et al, 2000). 3) In the close vicinity of the Galilean satellites (i.e. within several 
satellite radii) strong sharp peaks were seen in the dust impact rate. The spatial location 
of the events, as well as the impact velocities, establish a genetic link between the detected 
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dust grains and the satellites. Here we analyse this third population of dust in the close 
vicinity of the moons. First results have been published by Kriiger et al. (1999a). In this 
paper we present the detailed data analysis and the model for a dust cloud surrounding 
Ganymede. 

The most important question addressed here is: what is the physical mechanism that ejected 
the grains off the moons? Were these particles kicked up by hypervelocity impacts, pre- 
sumably of interplanetary micrometeoroids or tiny interstellar projectiles, onto the surface 
of the satellites? If this is true, we can treat the dust measurements made by Galileo at the 
Galilean satellites as a unique natural impact experiment which, for the first time, leads 
to the detection of the ejecta of hypervelocity impacts in space. This process has been 
suggested as being responsible for maintaining both the Jovian ring (Morfill et al., 1980; 
Hordnyi and Cravens, 1996; Ockert-Bell et al., 1999, Burns et al, 1999) and Saturn's E ring 
(Hordnyi et al, 1992; Hamilton and Burns, 1994). Our in-situ measurements in space give 
us a chance to gain more insight into the fundamental properties of hypervelocity impacts, 
for which the laboratory experiments on Earth still do not yield a comprehensive picture. 

Between December 1995 and April 2000 the Galileo spacecraft had a total of 26 targeted 
encounters with all 4 Galilean satellites: 4 with Io, 11 with Europa, 4 with Ganymede and 
7 with Callisto. In addition to the 4 Ganymede flybys which occurred in 1996 and 1997 
(encounters Gl, G2, G7 and G8) two more are planned for the year 2000 (G28 and G29). 
Due to the spacecraft orientation, however, the population of grains existing in the close 
vicinity of Ganymede will not be detectable during these two encounters. Hence, the data 
set of dust particles detected with Galileo during flybys at Ganymede will not be extended 
in the future. Furthermore, the analysis of the data obtained from the whole Galileo 
mission shows that the mass and velocity calibration of the dust instrument is subject to a 
radiation-related ageing effect which affects data obtained after 1996 and which is still under 
investigation. For these two reasons, we here focus our attention on the dust grains detected 
near Ganymede during the four close encounters with this satellite. Detailed investigations 
of the dust clouds surrounding Europa and Callisto will be the subject of a future paper. 

We begin with a description of the Galileo data and their processing (Section 2). Section 
3 presents the impact model that we use to interpret the data. A comparison of the data 
and modelling results is made in Section 4. Section 5 contains a summary and discussion 
of our findings. 



2 Dust impacts detected close to Ganymede 



Dust impact detection 

The dust detector onboard the Galileo spacecraft is a multi-coincidence impact ionisation 
dust detector which measures submicrometre- and micrometre-sized dust particles (Griin 
et al, 1992; 1995). The instrument is identical with the dust detector onboard Ulysses. 

For each dust grain hitting the detector target three independent measurements of the 
impact-created plasma cloud are used to derive the impact speed v and the mass m of the 
particle. The charge Qj released upon impact on the target depends on the mass and the 
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speed of the impacting grain according to: Qj oc m ■ i> 3 ' 5 . The calibrated velocity range 
of the instrument is 2 to 70kms _1 . The coincidence times of the three charge signals are 
used to classify each impact into one of four categories. Class 3 impacts have three charge 
signals, two are required for a class 2 and class 1 events, and only one for class 0. Class 3, 
our highest class, are real dust impacts and class are noise events. Class 1 and class 2 
events have been true dust impacts in interplanetary space (Baguhl et al., 1993, Kriiger et 
al, 1999c). 

Within about 15Rj distance from Jupiter energetic particles from the Jovian plasma envi- 
ronment cause an enhanced noise rate in class 2 and the lower quality classes. By analysing 
the properties of the Io stream particles and comparing them with the noise events, the 
noise could be eliminated from the class 2 data (Kriiger et al, 1999d), but class 1 shows 
signatures of being all noise in the Jovian environment. This algorithm, however, applies 
to the stream particles and cannot per se be applied to other populations of grains. Fur- 
thermore, the denoising technique uses statistical arguments and is applicable to large data 
sets. Individual dust impacts may be erroneously classified as noise and vice versa. The 
analysis of the data sets from the flybys of all three satellites (Europa, Ganymede and 
Callisto) implies a slightly modified algorithm for the rejection of noise events in the flyby 
data of Europa (Kriiger et al., 2000). 

To check for a possible contamination by noise in the data set used for our analysis, we 
take all class 3 (dust impacts) and all class 2 events (dust and noise) and compare this 
total number with the number of class 2 events classified as noise by the standard algorithm 
(Kriiger et al., 1999d) and by the modified algorithm applicable to the Europa data (Kriiger 
et al., 2000), respectively. A total of 27% of the events (10 out of 36) would be classified as 
noise in the first case and 14 % (5 out of 36) in the second one. Even the higher possible noise 
contamination in the first case does not affect the conclusions of this paper. In addition, 
Ganymede's orbital radius is at 15.0 Rj from Jupiter at the edge of the noisy region where 
the noise contamination is relatively small anyway. 

We therefore consider class 3 and all class 2 data for our analysis. It should be noted that 
there is no physical difference between dust impacts categorised into class 2 and class 3. 

Galileo is a dual-spinning spacecraft, with an antenna that points antiparallel to the space- 
craft positive spin axis. During most of the orbital tour around Jupiter, the antenna points 
towards Earth. The dust instrument is mounted on the spinning section of the spacecraft 
and its sensor axis is offset by an angle of 60° from the spin axis (Fig. |]; an angle of 55° 
was erroneously stated earlier). 

The rotation angle, 0, measures the viewing direction of the dust sensor at the time of 
particle impact. During one spin revolution of the spacecraft, the rotation angle scans 
through 360°. Rotation angles for the Galileo dust instrument, however, are reported 
opposite to that of the actual spacecraft rotation direction. This is done to easily compare 
Galileo results with the dust detector data taken on the Ulysses spacecraft, which, unlike 
Galileo, has the opposite spin direction. Zero degrees of rotation angle is taken when the 
dust sensor points close to the ecliptic north direction. At rotation angles of 90° and 270° 
the sensor axis lies nearly in the ecliptic plane (which corresponds roughly to Jupiter's 
equatorial plane) . The dust instrument has a 140° wide field of view (FOV) . Dust particles 
which arrive within 10° of the positive spin axis can be sensed at all rotation angles, while 
those that arrive at angles from 10° to 130° from the positive spin axis can only be sensed 
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over a limited range of rotation angles. 



Insert Figure [l] 

Because Galileo's high-gain antenna did not open completely, the spacecraft has a very 
low data transmission capability. For the dust instrument this means that the full set of 
parameters measured during a dust particle impact (rotation angle, impact charges, charge 
rise times, etc.) can only be transmitted to Earth for a limited number of particles. During 
the times considered in this paper these limits are between one particle per minute and one 
particle per 21 minutes (Kriiger et al., 2000). During periods of higher impact rates, the full 
set of parameters is not transmitted to Earth for all particles. All particles, however, are 
counted with one of 24 accumulators {Griin et at, 1995). This allows for the determination 
of reliable impact rates during the satellite flybys. 

Impact direction 

In Fig. |2] we show the impact direction (rotation angle) of those particles detected within 
about two hours around closest approach to Ganymede for which the complete set of mea- 
sured impact parameters has been transmitted to Earth. During the first three flybys (Gl, 
G2, G7) particles with rotation angles between 180° and 360° were strongly concentrated 
towards Ganymede. Most of them were detected at altitudes below 2 Rg (Ganymede radius 
Rq = 2,635 km). During the Gl and G2 flybys the large number of particles with rota- 
tion angles between 0° and 180° are streams of ten-nanometre-sized dust grains which have 
been detected throughout the Jovian system (see Griin et at, 1998, for a detailed analysis). 
The direction from which these dust streams were observed varied during Galileo's path 
through the Jovian system: when Galileo approached the inner Jovian system rotation 
angles around 270° were observed. Shortly before closest approach to Jupiter the rotation 
angle changed to 90°. Therefore, depending on when the satellite flyby occurred, stream 
particles approached from the corresponding direction: 90° in the case of Gl, G2 and G7, 
and 270° during G8. The stream particles can, in general, be identified by their calibrated 
impact velocities as will be shown further below. The detection geometry is sketched in 
Fig. |. 

At Gl and G2 an apparent concentration of particles with 0° < < 180° within about 
5 Rg altitude is due to a higher data transmission rate for about an hour around Ganymede 
closest approach (c/. Fig. ||). During this time period the number of particles for which 
their complete information has been transmitted to Earth is increased by about a factor of 
ten. At the G7 encounter almost all particles with 180° < < 360° were detected within 
4 Rg altitude. Only one stream particle was detected. During the G8 encounter the stream 
particles approached with rotation angles between 180° and 360°, and an apparent increase 
in the number of particles within about 4Rg altitude is again due to an enhanced data 
transmission rate. 

Insert Figure |^ 

The first three panels of Fig. ^ (Gl, G2 and G7) show that particles detected with 
180° < < 360° during the satellite flybys were concentrated towards Ganymede. They 
approached the dust sensor from a direction opposite to that of the dust stream particles. 
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We have analysed the velocity vector of Ganymede in a coordinate system fixed to Galileo 
in order to find out if the grains detected with 180° < < 360° are compatible with 
a Ganymede origin. This shows that particles approaching the sensor from Ganymede's 
direction which have a low velocity relative to Ganymede could be detected from a direction 
corresponding to about 270° rotation angle during all four encounters. Rotation angles of 
about 90° are not compatible with the Ganymede direction. Taking into account the sensor 
field of view of 140°, particles detected during the Gl, G2 and G7 encounters with rotation 
angles between 180° and 360° are compatible with an origin from Ganymede itself. In the 
following we will call them Ganymede particles. 

A total number of 36 Ganymede particles has been identified below 10 Rq altitude during 
these three encounters purely by their impact direction (c/. Tab. 1). For our further analysis 
we use a cut-off altitude of 10 Rq because this is close to the extension of the Hill sphere of 
Ganymede {Rh = 12 Rg)- F° r Gl and G2 the numbers of identified Ganymede particles are 
lower limits to the true numbers of detected grains because the complete set of parameters 
measured upon impact could be transmitted to Earth for only a fraction of all impacts. At 
G7, however, the complete set of parameters for all particles detected within a two hour 
period around closest approach has been transmitted. 

Insert Table 1 

Ganymede particles from the G8 encounter cannot be identified by their impact direction 
alone because both stream particles and Ganymede particles approached the sensor from 
the same direction. We therefore have to use an additional criterion to separate Ganymede 
particles at G8. 

Impact velocity 

If the Ganymede particles truly originate from the satellite itself and belong to a steady- 
state dust cloud they should impact the detector with roughly the velocity of the spacecraft 
relative to the moon. During all four flybys this encounter velocity was close to 8 km s _1 (see 
Tab. 1) which is well above the detection threshold of the dust instrument for micrometre- 
sized grains (2kms— 1). In contrast to the Ganymede particles, stream particles have 
velocities in excess of 200 kms -1 (Zook et al., 1996). Therefore, the impact velocity should 
be a good parameter to separate both populations of grains. 

Calibrated impact velocities are derived from the rise-time of the impact charge signal by an 
empirically derived algorithm. Such impact velocities could be determined for 29 Ganymede 
particles detected below 10 Rg altitude during the Gl, G2 and G7 flybys (velocity error 
factor VEF < 6; Griin et al., 1995). In Fig. || we compare the calibrated impact velocities 
of these Ganymede particles with those of the stream particles detected in the same time 
period. Because the distribution shows a bimodal structure, we can separate two statistical 
subsets: nearly all of the stream particles have calibrated velocities in excess of 10 kms -1 , 
whereas Ganymede particles are slower. The true impact velocities of the stream particles 
are much higher than the upper limit of the calibrated velocity range of the dust instrument 
(70 kms -1 ). Thus, the values derived from the instrument calibration are by far too low. 

Insert Figure ||| 
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The mean velocity of the 29 Ganymede particles in Fig. || is 7.2 ± 4.9 kms -1 (lcr). Given 
a typical uncertainty for an individual velocity measurement of a factor of two, this value 
is in good agreement with the velocity of Galileo relative to Ganymede (about 8 kms -1 ). 
It supports the picture that the particles detected belong to a dust cloud surrounding 
the moon. It also shows that the empirical velocity calibration of the instrument can be 
applied to the Ganymede particles, although one knows that it is not valid for the much 
smaller stream particles. Furthermore, and much more important for our analysis here, the 
calibrated impact velocities can be used to identify particles belonging to a dust cloud when 
such grains cannot be identified by their impact direction. It should be noted, however, 
that the two velocity distributions in Fig. |3| overlap which leads to some ambiguity in the 
identification of individual grains. 

In Fig. ^ we have marked particles according to their calibrated velocities: those with 
impact speeds below lOkms -1 are shown as circles, faster grains as crosses. For G8 (bottom 
panel) the majority of particles with 180° < G < 360° have calibrated impact velocities 
above 10 kms -1 : they are classified as stream particles and are rejected. Only 9 particles 
with 180° < < 360° remain below 5Rq altitude and a velocity below 10 kms -1 . Only 
particles with an altitude below 5Rq are considered for further analysis because we want 
to minimise the contamination by stream particles: a few stream particles in Fig. [3| have 
calibrated velocities below 10 kms -1 and would erroneously be classified as Ganymede 
particles. Such a cutoff altitude is implied by the absence of Ganymede particles during the 
G7 flyby when the full data set for all particles has been transmitted. The mean velocity of 
the 9 particles from G8 is 6.9 ± 2.4 kms -1 which is close to the value derived for the three 
earlier encounters. 

At the G8 encounter the identified Ganymede particles show a slight concentration towards 
the satellite. Together with the 36 particles identified by their impact direction during the 
three earlier flybys, we have identified a total set of 45 Ganymede particles from all four 
encounters. 

Impact rate 

Early analyses of the impact rate of dust particles measured with the Galileo dust detector 
during approaches to the Galilean moons showed a sharp peak within about half an hour 
of closest approach to the moon (cf. Griin et al., 1997, 1998; Kriiger et al., 1999b). This 
indicated the existence of concentrations of dust particles within a few satellite radii above 
the satellites' surfaces. 

With the data set for Ganymede particles detected during all four encounters we can now 
calculate the impact rate of dust grains in the close vicinity of the moon (Fig. . To obtain 
the impact rate we define distance bins equally spaced on a logarithmic scale and divide 
the number of particles for which the complete set of parameters has been transmitted to 
Earth in a given distance bin by the time Galileo has spent in this bin (dotted lines). 

Insert Figure ^ 

To correct for incomplete data transmission, we then multiply bin by bin the impact rate 
with the ratio between the number of counted particles and the number of particles for which 
the complete data set has been transmitted. The impact rates corrected for incomplete data 
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transmission are shown as solid lines. In cases where only a solid line is visible no correction 
was applied because the complete information of all particles in that bin is available (Gl 
and G7). 

During the first three encounters (Gl, G2 and G7) the impact rate increases towards 
Ganymede. Note specifically that for G7 no correction for incomplete transmission is needed 
because the full data set for all detected particles has been transmitted. Also, Gl needs 
only a correction in the innermost bins which steepens the slope slightly. Thus, Gl and G7 
probably give the most reliable slopes for the spatial distribution of the grains in the dust 
cloud. Both fits give power law slopes steeper than —2. 

For the fourth encounter (G8) no concentration of particles towards Ganymede is seen after 
correction for incomplete transmission. Note that a slight concentration is obvious in the 
uncorrected rate (dotted curve). This indicates that the slopes for the corrected curves have 
to be taken with some caution. Interestingly, the flattest slopes have been found for the 
encounters with the largest corrections for incomplete transmission (G2 and G8). On the 
other hand, the G2 flyby was the one which occurred closest to the north pole of Ganymede 
(80° latitude). We do not investigate variations of the slopes between individual encounters 
because of the large statistical uncertainties. Spatial variations with respect to the flyby 
position relative to the satellite will be addressed in a future investigation which will include 
data from Galileo's Europa and Callisto flybys. 

Mass distribution 

The charge released by an impact of a dust particle onto the target depends on the mass 
and the velocity of the grain (Griin et al., 1995). Thus, in order to calculate the particle 
mass from the measured impact parameters one has to know the impact velocity. In Fig. || 
we show the mass distribution of the particles from all four Ganymede flybys for which the 
velocity could be reliably determined (VEF < 6; 38 particles in total). 

Insert Figure |B| 

In the upper panel of Fig. ||| the impact velocities derived from the instrument calibration 
have been used to obtain the particle mass. With this method the uncertainty of the impact 
velocity is typically a factor of 2 and that of the mass is a factor of 10. 

The dust detector has a velocity-dependent detection threshold {Griin et al., 1995). The 
threshold for particles approaching with 8kms _1 is shown as a dashed line. The mass 
distribution is incomplete around and below this value. 

The mass distribution is also affected by the low data transmission capability of Galileo 
and the data storage scheme in the instrument memory. As a result, nearly all particles 
lost in Gl, G2 and G8 are in the mass range below 10 _15 kg. If we assume that the lost 
particles are equally distributed over the mass bins below this value, the maximum of the 
mass distribution is artificially too low by less than a factor of 2. 

If the individual impact velocities were known with a higher accuracy than the typical factor 
of 2 uncertainty from the instrument calibration, the uncertainty in the mass determination 
could be reduced. Because the mean values of the measured impact velocities are in good 
agreement with the velocities of Galileo relative to Ganymede for the individual encounters 
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(Tab. 1), we assume the latter ones as the particles' impact velocities and re-calculate the 
particle mass. This method has also successfully been applied to the mass distribution of 
interstellar dust particles (Landgraf et al., 2000). The lower panel of Fig. |H| shows the result 
for the Ganymede particles. The width of the mass distribution is significantly smaller than 
that derived from the calibrated impact velocities. 

The mass distribution of the grains allows for a simple check for the compatibility of the 
data with the hypothesis of the impact origin of the detected particles. We took the mass 
distributions shown in Fig. || and show linear fits to the cumulative distributions in Fig. ||. 
We only considered grains, the masses of which were definitely well above the detector 
threshold (M > 10 _16 kg). The slopes are in good agreement with the typical slopes one 
expects for the impact ejecta (0.5 & a 5s 1.0; see, e.g., Krivov and Jurewicz, 1999). 

Insert Figure ^ 



3 Model for the Impact Ejecta from Ganymede 

Detections of dust near Ganymede, and our results of the data processing which reveal 
the concentration of dust towards Ganymede, give rise to the following question: what 
is the physical mechanism that produces the dust material? A possibility of capturing 
the grains by the gravity of Ganymede is dynamically ruled out. In the 3-body problem 
"dust grain - Ganymede - Jupiter", a close encounter of a grain with Ganymede may 
convert a hyperbolic orbit with respect to Jupiter to an elliptical one, but still bound to 
Jupiter and not to Ganymede. Under some conditions, it is possible that interplanetary 
dust grains are trapped in Jovian magnetosphere by electromagnetic forces (Colwell et al, 
1998; Thiessenhusen et al., 2000); but a strong concentration of dust towards Ganymede 
requires additional mechanisms connected with Ganymede itself. The recently discovered 
own magnetic field of Ganymede {Kivelson et al., 1998) seems to be too weak to do the 
job. One could claim that geysers or volcanoes loft dust off Ganymede, but this is very 
unconvincing and unsupported by the observations which show no activity of this kind. 
One could argue for a self-sustained distribution of dust which was suggested to take place 
in Saturn's E ring {Hamilton and Burns, 1994), but the strong gravity of Ganymede makes 
that improbable. A conceivable guess one could make about the origin of the Ganymede 
grains is their production via continuous bombardment by interplanetary micrometeoroids. 

In what follows, we construct a model to predict the density of the dust cloud of Ganymede, 
produced by continuous hypervelocity impacts of interplanetary projectiles onto the moon's 
surface. We start with the main simplifying assumptions. Two populations of the debris 
particles are considered: those which move on ballistic trajectories and therefore fall back to 
the satellite shortly after they are ejected, as well as those which are fast enough to escape 
from Ganymede into circumjovian space. We assume that these ejecta move in Keplerian 
trajectories — pieces of ellipses and hyperbolas, respectively. This is justified by estimates 
of the perturbing forces. The tidal gravity of Jupiter is only important outside the Hill 
sphere (Rh ~ 12 Rq)- Next, for the grains with masses greater than 10~ 16 kg, to which we 
confine our analysis, the radiation pressure force is less than 10% of Ganymede's gravity 
even at 10 Rq from the moon. The Coriolis force is comparable to the satellite gravity at 
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~ 5Rq, but it cannot change appreciably the radial distribution of the ejecta. The only 
force that can possibly make the motion essentially non-Keplerian is the Lorentz force. Still 
this force, which is very important for 10~ 16 kg grains moving about Jupiter (Colwell et al., 
1998; Krivov et al., 2000), can probably be neglected for short-living grains in the vicinity of 
Ganymede. Assuming a strength of the Jovian magnetic field at Ganymede of ~ 10 2 nT, the 
velocity of the grains to be equal to the Ganymede velocity relative to the corotating jovian 
magnetic field (180 km s -1 ), and electrostatic potentials of grains of +3V (Hordnyi et al., 
1993), we find that at 5Rg the Lorentz force amounts to ~ 30% of the satellite gravity. The 
intrinsic magnetic field of Ganymede is ten times stronger (£> 10 3 nT; Kivelson et al, 1998), 
but the velocity of the grains relative to Ganymede and its magnetic field is only of the 
order of 1 kms -1 , which means that the relative strength of the Lorentz force caused by the 
satellite's magnetic field is one order of magnitude weaker. We conclude, especially taking 
into account the uncertainty in grains' charges, that under some conditions the Lorentz 
force might lead to noticeable effects. However, incorporation of the electromagnetic effects 
into the model would require an extensive dedicated modelling, including implementation 
of models of Ganymede's magnetosphere and performing the charging calculations. This 
would make little sense in view of the scarcity of the Galileo dust data. Thus this task is 
beyond the scope of the present paper. 

We also assume that the circumsatellite dust cloud is in a steady state. Finally, we neglect 
possible effects of non-isotropy of the impactor flux (e.g., Colwell, 1993) and assume that 
the cloud is spherically symmetric. Thus, in this paper we do not consider variations in the 
spatial density from flyby to flyby which may be caused by spatial or temporal variations 
of the dust cloud surrounding Ganymede. Again, such an investigation is not possible with 
the present data set of 45 ejecta grains. 

Impactor flux and ejecta yield 

To predict the flux of the interplanetary impactors, we use Divine's (1993) model. We have 
calculated the fluxes of interplanetary grains onto a sphere with unit cross section (irR 2 = 1), 
moving around the Sun in a circular Keplerian orbit with a radius of 5.2 AU (heliocentric 
distance of Jupiter). Figure |7| shows, as a function of particle's mass, the cumulative flux 
onto this sphere, the differential flux per mass decade, as well as the differential mass flux 
per mass decade. For comparison, the same fluxes at 1AU are also shown. The total 
mass flux at 5.2AU "far" from Jupiter (i.e., before its gravity is taken into account) is 
Foo = 8.4 x 10~ 16 kgm 2 s . It is dominated by grains with masses ~ 10 -8 kg (or with sizes 
~ 100 (tm). 

Insert Figure |7] 

The main contribution to the flux and mass flux is made by the "halo population" (Di- 
vine, 1993), composed of particles in orbits with moderate eccentricities, but randomly 
distributed inclinations between 0° and 180° (thus including grains in both prograde and 
retrograde orbits). The velocities of interplanetary grains with respect to the Sun are very 
high: ~ 20 kms . Their mean velocity relative to Jupiter, because of the random distri- 
bution of inclinations, is nearly the same, ~ 20 kms -1 . Obviously, the dispersion is also 
high: since the orbital velocity of Jupiter round the Sun is ~ 13 kms -1 , the velocities of 
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individual particles may range from 7 to 33kms _1 . For comparison, Colwell and Hordnyi 
(1996) adopt for this "Oort cloud population" of grains the velocity range of 8 to 34kms _1 
with a mean of 23.6 km s -1 . 

We have to take into account the gravitational focussing by Jupiter, which increases the 
interplanetary flux F at Ganymede (r = 15Rj): (i) the speed of the grains v becomes larger 
than that far from the planet, Voo, and (ii) the spatial density of dust n gets larger than the 
one far from Jupiter, tIqq. Assuming = 20 km s 1 and applying the energy integral, we 
find v/voo = 1-3 and v ~ 25kms~ 1 . The velocity v has the meaning of the mean velocity of 
impactors with respect to Jupiter at Ganymede's distance. Nevertheless, since the incoming 
directions of the projectiles are broadly distributed, we can roughly take this value as the 
mean projectile velocity with respect to Ganymede. We note, however, that the velocities of 
individual grains striking Ganymede may range from nearly zero to about 50kms _1 . The 
former takes place when Ganymede is at opposition to the Sun for the prograde particles 
hitting its trailing hemisphere; the latter is attained at the same moment for retrograde 
particles impacting the leading hemisphere. Next, using the formulae by Colombo et al. 
(1966), we find nju^ = 1.4. Therefore, 

F v n 

-ET- = = 1-8, (1) 

so that the mass flux at Ganymede is F = 1.5 x 10 -15 kgm 2 s _1 . 

We consider the impact ejecta production now. The efficiency of the material ejection in 
a cratering event is characterised by Y, the characteristic yield, defined as the ratio of the 
ejected mass to the projectile mass. For hypervelocity impacts into ices, the typical yields 
Y were reported to range from ~ 10 3 to 10 6 (e.g., Lange and Ahrens, 1987; Koschny, priv. 
comm.) .From (Koschny, priv. comm.), for pure ice, for 10~ 8 kg impactors, and for an impact 
speed of 25kms _1 , the yield is Y ~ 1 x 10 4 . The mass production rate from Ganymede's 
surface is then 

M + = FYS = 330kgs- 1 , (2) 

where S = irR^ = 2.2 x 10 13 m 2 is the cross section area of Ganymede. 

The mass distribution of the ejecta is commonly represented by a power law: iV + (> M) oc 
M~ a , where M is the grain mass and N + (> M) is the number of particles with masses 
> M ejected from the moon per second. The distribution is normalised to the total amount 
of the ejecta M+: 

1-a M+ (M max 



AT+(>M)= ^ , (3) 

a M max \ M J 

where M + is the largest mass of the ejecta. A plausible slope of the ejecta mass distribution 
is a ~ 0.83 (e.g., Koschny, priv. comm.; Krivov and Jurewicz, 1999) which is consistent with 
the exponent found from the Galileo data (Fig. ||). For the heaviest ejecta fragment, we take 
M-max = 10~ 8 kg, which is close to the typical mass of the impactors (the dependence of this 
parameter is weak though). In what follows, we always calculate = iV + (> 10~ 16 kg), 
which correspond to the grains "safely" above the detection threshold, cf. Fig. ||). 



Ejecta speed distribution 
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In this study, we assume no dependence between the speeds of the ejected grains and their 
masses, because laboratory experiments still do not evidence any strong correlation between 
the two quantities (e.g., Nakamura and Fujiwara, 1991). Thus, regardless of how massive 
the ejecta are, we take the ejecta speed distribution in the form (see, e.g., Stoffler et al., 
1975, Hartmann, 1985) 

*(> u) = (u/u )-\ (4) 
where ^(> n) is the fraction of the material ejected at speeds > n. 

This adds two parameters to the model: the "lower cut-off velocity" no and the distribution 
slope 7. The values of these parameters were determined in impact experiments. A plausible 
range for no could be from a few ms" 1 to several hundred ms _1 (e.g., Stoffler et al., 1975; 
Hartmann, 1985). The slope of the speed distribution 7 may range from 1 to somewhat 
greater than 2 (e.g., Frisch, 1992; see also a discussion in Colwell and Esposito, 1993). 

A remark should be made about the directions of the ejecta velocities. From the labo- 
ratory experiments it is known that even oblique impacts lead to ejection of the surface 
material into a cone normal to a target surface (Nakamura and Fujiwara, 1991). By nu- 
merical experiments, we have checked that the spatial distribution of the ejecta depends 
only weakly on the opening angle of this cone. Also, laboratory measurements suggest 
that most of the debris are ejected at high angles to the surface — such as 60° (Koschny, 
priv. comm.;Burchell et al., 1998; also M. J. Burchell, priv. comm.) or even 80° (Frisch, 
1992). Besides, some authors observed vertical plumes of fast ejecta (Hartmann, 1985). For 
these reasons, and to alleviate the analytical derivations, we assume for simplicity that the 
trajectories of the ejecta are straight lines normal to the surface of Ganymede ("degenerate 
Keplerian ellipses"). 

Energy constraints on the model parameters 

The model as a whole contains five parameters. Three parameters come from the ejecta 
mass distribution (characteristic yield Y, slope a, and the maximum mass of a fragment 
M max ), while two others pertain to the ejecta speed distribution (the lower cut-off velocity 
no and the slope 7). The dependence of the results on a and M max is only moderate, and 
besides, changes in their values can always be absorbed in a value of the characteristic yield. 
This leaves Y, no, and 7 as the three most important model parameters. Of course, the 
parameters describing the projectile flux, most notably the mass flux of the impactors F, 
are implicit in the results. 

Some constraints on the parameters Y, uq, and 7 can be placed by considering the energy 
balance. The ratio of the kinetic energy carried by the ejecta cloud, K e , to the impactor's 
kinetic energy, Ki, must be less than unity, because part of Ki is spent for comminution 
and heating. The ratios determined in the hypervelocity impact experiments vary from 
~ 0.2 to ~ 0.5 (e.g., Asada, 1985; Hartmann, 1985). We assume / = K e /Ki = 0.3. 

Let us express Ki and K e through the parameters of the model. Obviously, Ki = mv 2 /2, 
and 

K e = Ym —^(u)du, (5) 

where m is the impactor mass, v is its speed, u max is the speed of the fastest ejecta 
fragments, and V( M ) = — ^'(> n) (hereafter prime means derivative) is the differential 
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speed distribution, with ty{> u) being the cumulative distribution given by Eq. (||). We 
have: 

7-2 

- 1 



f = Y 



7 (Uq 



2 — 7 V v 



"o 



(7 ¥= 2). 



(6) 



(If 7 = 2, a logarithmic dependence results.) In Fig. ||, the lines / = 0.3 are plotted on the 
7, uo-plane for several values of Y. For Y = 10 4 , which we have chosen as the most plausible 
value, and for 7 = 1.2, 1.6, and 2.0, one gets uq = 13, 30, and 48ms , respectively. In 
these calculations, we assumed u max = 3kms _1 . The results do not depend on this speed 
severely. For example, with u max = 4kms _1 , for the same Y and 7 one gets uq = 10, 27, 
and 45 ms -1 . The values of uq that we derived are consistent with impact experiment data 
{e.g., Stoffler et al, 1975; Hartmann, 1985). 



Insert Figure 



Number density of dust 

Given a dust production rate N + , which can be calculated using Eq. (||), and assuming the 
ejecta speed distribution (Eq. ^[), we compute now the steady-state number density of the 
ejecta, n, as a function of r, the distance from the centre of Ganymede. 

Consider first the ejecta in bound orbits, i.e., the grains ejected at speeds u < u esc , where 
u esc = 2,750 ms" 1 is the escape velocity from Ganymede's surface. Remember that we 
assume the grains to move in rectilinear Keplerian trajectories. As follows from Eq. (0) and 
from the energy integral of the two-body problem, the fraction of ejecta, ejected in orbits 
with semimajor axes > a, is 

*(> a) = («o/««c) 7 (l - R G /{2a))-</ 2 . (7) 

The steady-state number of grains with semimajor axes (a, a + da) at distances (r, r + dr) 
from Ganymede's centre, where r < 2a, equals 

- 2iV+^'(> a)dadt, (8) 

where the factor 2 appears because a grain reaches the distance r twice — moving from and 
back to the moon, and dt is the time interval it takes for a grain to move from r to r + dr. 
Dividing Eq. (||) by the volume 4nr 2 dr of a spherical layer (r, r + dr) and integrating over 
semimajor axes of grains that reach distance r, one finds the steady-state number density 
of grains at distance r: 

j\f+ r+00 fa 
n bound {r) = ^ J (-*'(> a))-da. (9) 

From the energy integral, 

* = J_ / r l R o m 
dr u esc \ l-r/(2a)' V ; 
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Inserting Eq. (Q) and Eq. (|i~0|) into Eq. and using a new integration variable x = 2a/r, 
after some algebra we find 



-oc 



_ N+_ _j_ fuo_y Rg f dx 

b ° Und[ ) ~ ^ u esc \u esc ) V r 'J (l-^r-iz-i^xVT^T- ' ' 

Now we perform a similar derivation for the escaping grains with velocities u ^ t^esc 

at 

the surface. The steady-state number of grains ejected at speeds (u, u + du) (u > u esc ) at 
distances (r, r + dr) from Ganymede's centre is 

- N + V'(> u)dudt, (12) 



where dt is again the time that a grain needs to move from r to r + dr. Dividing Eq. (|1 
by the volume of a spherical layer (r, r + dr) and integrating over u > u esc there results: 

]\f+ r+OO fif 
n U nbound(r) = J— % J (~* (> u)) — du. ( 13 ) 

Here, 

* = ^ (14) 

dr u esc yj(u/u esc ) 2 - 1 + .Rc/r 

Substituting Eq. (|) and Eq. (0) into Eq. |0J), and using an integration variable x — ufu esc ^ 
after some transformations one obtains 

AT+ 7 /uoV T dx 
n unb ound(r) = — o- — • / ^— , .. ===== . ■ (15) 



4vrr 2 u esc \u esc / 7 x i+ 7 ^/ x 2 _ J + RG j 



r 



Let us discuss some properties of the functions nb OU nd( r ) (Eq. [Tl]) and n un i, oun( i(r) (Eq. |i~5| ). 
To get the actual number density at a distance r, we sum up the two contributions: n(r) = 
nbound{r) + n unbound (r). As seen from Eqs. (||) and (|l5|), n bound (r) decreases faster with 
increasing r than n unbound {r) , so that the relative contribution of escaping grains increases 
with distance. However, how large is this contribution at a given distance from the satellite? 
For one set of the model parameters, namely Y = 1 x 10 4 and uq = 30 ms -1 , we made 



calculations by using both Eqs. (11) and (15) and found that the number densities of 



the escaping grains and the particles falling back to the moon become comparable at the 
altitude of ~ 8Rq from Ganymede (Fig. ^a). As the Galileo detections have been made 
closer to the moon, we can conclude that the majority of the ejecta particles identified by 
the dust instrument, if they had not been caught by the detector, eventually would have 
fallen back onto Ganymede. 

Somewhat counterintuitively, the radial slope of the number density is not very sensitive to 
7. For all reasonable values of the ejecta distribution slope 7, the number density is close 
to a power law n(r) oc r~~ u with v ~ 2 to 3 (see Fig. ||b for grains in bound orbits and 
Fig. ||c for escaping grains). 

Remember that, in derivation of Eqs. (|ll|) and ([15])) we assumed the ejecta to move radially 
from and to the surface of Ganymede. This assumption is partly reflected by the lower 
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integration limit of 1 in Eq. (|ll]). We calculated the integral with integration limits greater 
than unity, to simulate the particles leaving the surface not vertically. The results (Fig. ^d) 
show that at small altitudes the number density gets smaller as compared to the case of 
vertical ejections, and the curves slope more gently. 

Insert Figure ^ 



4 Comparison of the Model with Galileo Data 

Comparison of the available data with theory is accomplished by constructing the depen- 
dencies of the frequency of the detections as a function of radial distance from the moon. 
To reach this goal, we use the impact rates depicted by solid lines in Fig. |j. Dividing the 
rates by the effective spin-averaged detector area, we obtain the fluxes (m _2 s _1 ). Then 
we divide the results by the mean impact velocity for a given flyby, which results in mean 
number densities (m -3 ) in various distance bins explained in Sect. 2. 

On the other hand, we use our model to predict the number density as a function of distance. 
We used Eqs. (pi), (11) and (15). As discussed above, we chose the yield Y = 10 4 and took 



parameters no, and 7 in the ranges compatible with laboratory impact experiments, using 
the values that satisfy the relation K e /Ki = 0.3. As regards two other, less important, 
model parameters, the values a = 0.83 and M max = 10~ 8 kg (cf. Asada, 1985; O'Keefe and 
Ahrens, 1985; Koschny, priv. comm.) were used. 



The results (both the data points and theoretical curves) are displayed in Fig. 1C, showing 
a reasonably good agreement between the data and the model. For example, Y = 10 4 , 
no £ [13,48] ms^ 1 , and 7 E [1.2,2.0] are quite compatible with the data for Gl and G7, 
which represent the most reliable number density slopes. We cannot leave without notice, 
however, that the theoretical curves seem to be somewhat steeper than the distribution of 
the data points for G2 and G8. Either this effect is caused by statistical errors (see Section 
2), or it may be attributed to some simplifying assumptions of our model (see Sections 
3 and 5). Note that using the lower integration limit in Eq. (|i~T| ) larger than 1, which 
simulates the oblique ejections from the surface, decreases the predicted number density 
at lower altitudes, resulting in a better compatibility with some of the data points (e.g., 
leftmost G2 point in Fig. 10). 



Insert Figure 



It should be emphasised that in the present paper, we do not make any systematic attempt 
to constrain the poorly known parameters of the model from the data. One reason for that 
is scarcity of the data. Another one is that some of the parameters (such as a, M max , 
and 7) affect the modelling results only slightly, while the others (such as Y, no, and 
the parameters describing the impactor flux) are strongly cross-correlated and cannot be 
constrained independently. A weak dependence of the predicted number density on some 
of the parameters would certainly be a serious handicap, if our goal were to retrieve the 
poorly known parameters. However, the same circumstance — a weak sensitivity of the 
model results to the model parameters — turns to a serious advantage, if we are merely 
checking the compatibility of the data with our model. The fact that the model agrees 
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with the data over a wide range of the parameters strengthens our conclusion about the 
likelihood of the impact origin of the detected grains. 

In addition to the study of the debris produced by impacts of interplanetary grains, we have 
also estimated a possible effect of interstellar grains as impactors. The result is that their 
contribution is negligible. The reason is that, although the flux of interstellar grains at the 
heliocentric distance of Jupiter is known to exceed that of the interplanetary particles, the 
mass flux in the considered mass range is about 4 orders of magnitude smaller (Landgraf et 
al, 2000). Hence the flux of the ejecta produced by interstellar grains is also substantially 
smaller than the one caused by interplanetary projectiles. 

5 Conclusions and Discussion 

In this paper we examined the dust impacts registered by the Galileo dust detector in the 
immediate vicinity of Ganymede during four close flybys of this Jovian moon. Analyzing 
impact directions and velocities and the mass distribution, as well as spatial locations of the 
events, we have evidenced that the particles did originate from the moon. In an attempt 
to find out the specific physical mechanism that ejected the particles off the surface, we 
checked whether these "Ganymede dust grains" are impact debris produced by hyperve- 
locity impacts and, which impactors are responsible for the ejections. To do so, we have 
constructed a model for the impact ejecta from Ganymede and compared the modelling re- 
sults with the data. The comparison made it clear that the data are fully compatible with 
the impact origin of the Ganymede dust grains. Of course, this result cannot be treated 
as an evidence of the impact origin of the detected dust. With our modelling, we have 
only demonstrated that the data do not contravene the impact scenario, given our current 
knowledge of the physical conditions in the Jovian system as well as available laboratory 
data of hypervelocity impacts. 

If the impact scenario is true, then the class of projectiles responsible for the formation of 
a circumsatellite cloud of dust grains, several tens of which were sensed by the detector 
onboard Galileo, is almost definitely interplanetary micrometeoroids. The derived mass 
distribution of the detected grains, as well as the impact rates actually measured by the dust 
instrument, are in fairly good agreement with what we get from the model of hypervelocity 
impacts of interplanetary dust particles (IDPs), assuming contemporary models of IDP 
flux at a heliocentric distance of Jupiter and a low-temperature ice target. The interstellar 
grains should be much less efficient as impactors in producing the collisional debris, because 
of the much lower mass influx of interstellar dust onto the Ganymede surface compared to 
that of the IDPs. 

As follows from elementary estimates, the particles that reach an altitude of one satellite 
radius must have a starting velocity in excess of 2kms _1 . This confirms laboratory results 
that some of the hypervelocity impact ejecta from icy targets attain very high speeds 
(Frisch, 1992). Nevertheless, more than a half of the grains which could be found (and have 
been actually detected by Galileo) at altitudes less than about eight satellite radii above the 
surface, are not the particles escaping into Jovian space. They are slower ejecta, destined 
to fall back to Ganymede typically within several minutes to several hours after ejection. 
Such grains surround the satellite all the time as a result of continuous bombardment of 
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the surface by IDPs. Further out, starting from distances of ~ 8Rg, the escaping grains 
start to dominate the number density. 

It seems useful to give some general estimates concerning the mass budget of the dust cloud 
of Ganymede. The mass flux of IDPs bombarding the satellite surface is estimated to be 
~ 30 gs -1 (dominated by IDPs with m ~ 1CT 8 kg). With the characteristic yields of ~ 10 4 , 
we then estimate that as much as ~ 10 2 to lC^kgs- 1 of the moon's surface material, is 
ejected into space. Depending on the ejecta speed distribution adopted, the mean lifetime 
of the ejecta ranges from tens of seconds to several hundred seconds. Note that these values 
are dominated by the slowest ejecta — for the grains that reach the altitudes of several 
satellite radii the flight times are rather several hours. We therefore get an estimate of the 
total amount of dust contained in a steady-state cloud around Ganymede: 1 to 100 tons. 
Next, about 10~ 6 to 10~ 3 of the ejecta escape from the satellite. It gives the injection rate 
of the material into circumjovian space of ~ 10~ 4 to lkgs -1 , which may be comparable 
with the influx rate of IDPs to the Ganymede surface. 

Remember that, in our data interpretation, we assumed no time variations of the incoming 
impactor flux, so that the dust cloud was assumed to be in a steady state. It would be 
interesting to lift this assumption and to estimate possible time variations — e.g., tempo- 
rary enhancements of number density caused by hits of more massive meteorites, treated 
as individual events. It should be noted, however, that no pronounced enhancements of 
concentration of micrometre- sized debris are expected (Krivov and Jurewicz, 1999). 

Another, and possibly more restrictive simplification that we made is the spherical symme- 
try of the dust cloud, which implies an isotropy of the impactor flux. The mean speed of 
IDPs with respect to Jupiter at the distance of Ganymede is about 25 km s _1 (see Section 3). 
Jupiter's motion about the Sun (with the velocity of 13kms~ 4 ) and the orbital motion of 
the moon about Jupiter (11 km s -1 ) make the distribution of impactors at Ganymede highly 
unisotropic. Considering the geometry of the orbital motions of Jupiter, Ganymede, and 
IDPs, it is easy to show that velocities of the grains with respect to Ganymede vary from 
nearly zero to about SOkms^" 1 , depending on the position of Ganymede on its orbit and the 
directions of the particle velocities. As explained in Section 3, these extreme velocity values 
are attained when Ganymede is at opposition to the Sun, at which moment we can expect 
a very high asymmetry of the dust cloud above the leading and trailing hemispheres of 
the moon. The effect gets smaller at Ganymede's conjunction, when the prograde particles 
hit the trailing hemisphere with > 23 and retrograde grains reach the leading side with 
< 27kms~ 1 , but still many more impacts are expected on the leading edge of Ganymede 
than on the trailing edge. Therefore, we predict a marked and time-variable hemispheric 
asymmetry. Unfortunately, given the large statistical uncertainties due to the small num- 
ber of detected particles, no obvious spatial and/or temporal variation could be found in 
the present data. For this reason, we do not undertake detailed modelling for the angular 
distribution of impactors and ejecta similar to what has been done in different problems 
(e.g., Zook, 1992; Colwell, 1993). We hope to address these issues in future investigations 
of the full data set for all Galilean moons. 

Much of the dust ejected from Ganymede is launched into bound orbits and falls back to the 
moon. A complex of these short-living, but continuously replenished grains forms what we 
call an ejecta dust cloud of Ganymede. Obviously, all massive satellites which lack gaseous 
atmospheres should own an ejecta dust cloud. Before Galileo, there were few attempts of 
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direct in-situ detections of large ejecta close to satellites — most notably, near the Moon 
(Iglseder et al., 1996); these experiments did not bring up definite results, however. Here 
we present the successful measurements of the satellite ejecta in the vicinity of a source 
moon. We also present a dynamical model for an ejecta dust cloud of a massive satellite. 

A tiny fraction of impact debris is ejected at speeds sufficient to escape from Ganymede 
entirely. As shown above - within the error bars - the ejected mass is comparable with the 
incoming flux of IDP imactors. The ejected material goes into orbit around Jupiter and 
most of it will eventually be swept up by one of the Galilean satellites. These escaping grains 
are probably responsible for some of the impact events detected by the dust instrument in 
the inner Jovian system between the Galilean satellites (Griin et al, 1997, 1998; Krivov 
et al. 2000). Unfortunately, the ring of material formed by these grains escaping from 
Ganymede is far too tenuous to be detected optically. However, the fraction of debris 
escaping a satellite is a steeply decreasing function of satellite mass, so steep that despite 
their reduced cross sections, small moons may be better sources of dust than large satellites. 
Indeed, many small moons have been proved to be, or are supposed to be, sources of material 
for tenuous dusty rings surrounding all giant planets (Burns et al., 1984) and presumably 
Mars (Soter, 1971; Krivov and Hamilton, 1997). In the Jovian system, four innermost 
moons are sources for the dust halo, the main ring, and the two-component gossamer ring 
(Ockert-Bell et al., 1999; Burns et al., 1999). A variety of examples can be found in the 
Saturnian system, which is especially important in view of the Cassini mission. Enceladus 
supplies material to the huge E ring (Hamilton and Burns, 1994). Irregularly-shaped icy 
Hyperion is believed to be a source of icy ejecta which arrive at Titan, possibly affecting 
the complex chemistry of its dense nitrogen atmosphere (Banaszkiewicz and Krivov, 1997; 
Krivov and Banaszkiewicz, 2000). Yet farther out from Saturn, the outermost retrograde 
moon Phoebe emits dust which is thought to be deposited on the leading side of Iapetus, 
producing its observed brightness asymmetry (Burns et al., 1996). All these cases exemplify 
the same mechanism of dust production as the one discussed here for Ganymede. 

Therefore, spacecraft measurements near the satellites, i.e. very close to the sources of 
dust, would be of primary importance to gain more insight into the properties of satellite 
surfaces and the dusty rings these moons maintain. We believe the technique of processing 
and interpretation of spacecraft data which we proposed here and tested on Ganymede can 
be directly applied to data from the dust experiment onboard the Cassini spacecraft which 
are to be obtained during its fiybys of the Saturnian moons. 

The Galileo dust measurements at Ganymede presented in this paper can be considered 
as a unique natural impact experiment. They complement laboratory experiments in an 
astrophysically relevant environment. Laboratory impact experiments have significant de- 
ficiencies in many respects, in the speeds of the projectiles and the mass and speed ranges 
in which ejecta particles can be observed. Furthermore, there is the always pending ques- 
tion of the astrophysical relevance of the materials used. Although far from being perfect 
impact experiments, the Galileo results offer two extremely important improvements over 
laboratory experiments: 1) the projectile and target materials and projectile speeds are 
astrophysically relevant, and 2) the masses and speeds of the ejecta particles can be deter- 
mined in an important region of parameter space (micrometre sizes and kms -1 speeds). 
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Table 1: Parameters for the dust particles detected within 10 Rq altitude during Galileo's 
Ganymede flybys. The first two columns list the flyby number and the time of the flyby 
(year-day of year), column 3 gives the altitude at the closest approach to Ganymede, column 
4 gives the velocity of Galileo relative to Ganymede, column 5 lists the average particle 
velocity (only particles with a velocity error factor VEF < 6 have been used, cf. Griin et 
al. (1995), column 6 gives the number of class 2 and class 3 Ganymede particles for which 
the complete data set has been transmitted to Earth, and column 7 gives the number of 
Ganymede particles after correction for incomplete data transmission and column 8 the 
slope of the power law fit to the radial variation of the impact rate (cf. Fig. particles 
below 10 Rq altitude have been included in the power law fits for Gl, G2 and G7, and 
below 6Rq for G8, respectively). 



Flyby 


Date 


Altitude 


Spacecraft 


Average 


Particles 


Corrected 


Slope of 




(year-day) 




velocity 


particle 


with full 


number of 


impact 










velocity 


data set 


particles 


rate 






(km) 


( kms" 1 ) 


(kms -1 ) 








(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


Gl 


96-179.270 


844 


7.8 


6.3 ±2.1 


15 


30 


-2.45 ±0.16 


G2 


96-250.791 


262 


8.0 


8.7 ±7.8 


10 


48 


-0.52 ±0.10 


G7 


97-095.299 


3,095 


8.5 


6.7 ±3.7 


11 


11 


-2.32 ±0.02 


G8 


97-127.665 


1,596 


8.6 


6.9 ±2.4 


9* 


49 


0.37 ±0.22 



*Only particles with impact velocity v < 10 kms 1 and below 6Rg altitude included. 
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Figure 1: Galileo's trajectory and geometry of dust detection during the Gl Ganymede 
flyby. The Galileo spacecraft is sketched in an orientation it was in during the flyby (see 
text for details). The directions to Jupiter, Earth and Sun are shown. C/A indicates closest 
approach to Ganymede, FOV the field of view of the dust instrument. The orientation of 
the dust instrument shown corresponds to a rotation angle 6 = 270° . At 90° rotation angle 
it points in the opposite direction. 
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Figure 2: Sensor direction (rotation angle, O) versus altitude of the Galileo spacecraft 
above the surface of Ganymede at the time of dust impact. Data are shown for all four 
Ganymede encounters (Gl, G2, G7, G8). The radius of Ganymede is Rq = 2,635 km. The 
altitude range shown corresponds to a time interval of 2 h. Each symbol indicates a dust 
particle impact and the size of the circle indicates the impact charge created by the particle 
(10~ 14 C < Qi < 10~ n C). Circles show particles with impact velocities below lOkms -1 
and crosses show particles with higher speeds. Galileo did not traverse the region between 
the vertical dashed lines. 
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Figure 3: Impact velocities derived from the instrument calibration for dust particles 
detected at the Gl, G2 and G7 encounters below 10 Rg altitude. The solid line shows the 
distribution for Ganymede particles (180° < < 360°) and the dotted line that for stream 
particles (0° < < 180°). The mean velocity is 7.2 ±4.9 kms -1 . Only particles with 
a velocity error factor VEF < 6 (Griin et al, 1995) have been considered (29 Ganymede 
particles) . 
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Figure 4: Impact rate of dust particles detected during the four Ganymede encounters. 
The dotted histogram bins show the impact rate derived from the number of particles for 
which their complete information has been transmitted to Earth. The solid histograms 
show the same rates, but corrected for incomplete data transmission. The vertical dotted 
lines indicate the minimum altitude reached by Galileo at closest approach. Error bars 
denote the y/n statistical uncertainty, with n being the number of particles for which the 
complete information has been transmitted. For G8 only particles with a calibrated velocity 
below 10 kms" 1 have been used, and the rate has been truncated at 6Rq- Dashed lines 
are power law fits to the corrected impact rate (c/. Tab. 1). 
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Figure 5: Mass distribution of the particles detected during the four Ganymede flybys. 
The upper panel shows the distribution obtained by using the measured impact velocities 
derived from the instrument calibration. In the lower panel the velocity of Galileo relative 
to Ganymede has been assumed as the impact velocity in order to calculate the mass of the 
particle. The dashed lines indicate the detection threshold for particles which approach the 
detector with the velocity of Galileo relative to Ganymede (about 8kms _1 ). Only the 38 
particles with a velocity error factor VEF < 6 (Griin et ai, 1995) have been considered. 
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Figure 6: Cumulative mass distributions of Fig. ^ and linear fits to the data. Vertical bars 
indicate the \fn statistical uncertainty. 
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Figure 7: Fluxes of interplanetary grains onto a sphere with unit cross section (ttR 2 = 1), 
moving around the Sun in a circular Keplerian orbit with radius of 5.2 AU (solid lines) and 
1 AU (dashed lines), according to Divine's (1993) model. Top: the cumulative flux; middle: 
the differential flux per unit log mass interval; bottom: the differential mass flux per unit 
log mass interval. 
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Figure 8: Constraints on the model parameters Y, uq, and 7 coming from the consideration 
of the kinetic energy. Shown, for several values of the characteristic yield Y, are the lines 
of constant ratio K e /K{ = 0.3, where Ki and K e are respectively the kinetic energy of 
an impactor and the ejecta it produces. The impact speed is taken to be 25kms _1 . The 
speeds of the ejecta are assumed to obey a power law (Eq. ^) with a slope 7, a lower cut-off 
uq and the upper cut-off u max = 3kms _1 . 
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Figure 9: Plots of the functions nf, ou „d(r) (Eq. [Tl]) and n un b oun( i(r) (Eq. |15|). Adopted 
values of the parameters: Y = 1 X 10 4 , uq = 30 ms -1 . (a) These two functions and their 
sum for a fixed 7 = 1.6. (b) Number density of grains in bound orbits for several values 
of 7. (c) The same, but for escaping grains, (d) Number density of dust (in bound and 
unbound orbits together) for a fixed 7 = 1.6, but for different lower integration limits in 
Eq. (11), which simulates oblique ejections of grains from the surface. See text for further 
explanations. 
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Figure 10: The number density of dust as a function of distance from the centre of 
Ganymede, derived from the data (symbols with error bars) and predicted by the model 
(lines). Horizontal bars for the data symbols indicate distance bins which were used in 
processing the data (see text for details), whereas vertical ones reflect ^/n errors due to 
a limited number of impacts. Theoretical curves are given for several plausible choices of 
the model parameters, listed in the legend. For one of the models, we show the result 
for simulated non- vertical ejections of dust from the surface (the lower limit in Eq. ( |ll| ) 
was taken to be 1.25 instead of 1). All the lines show the total number density of ejecta 
(falling back and escaping) produced by interplanetary impactors. Here we do not consider 
variations in the spatial density from fiyby to flyby which may be caused by spatial or 
temporal variations of the dust cloud surrounding Ganymede. 
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